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Abstract 

The finite-volume nature of lattice QCD entails a variety of effects that must be 
handled in the process of performing chiral extrapolations. Since the pion cloud 
that surrounds hadrons becomes distorted in a finite volume, hadronic observables 
must be corrected before one can compare with the experimental values. The 
electric charge radius of the nucleon is of particular interest when considering 
the implementation of finite-volume corrections. It is common practice in the 
p I literature to transform electric form factors from the lattice into charge radii prior 

O . to analysis. However, there is a fundamental difficulty with using these charge 

radii in a finite-volume extrapolation. The subtleties are a consequence of the 
absence of a continuous derivative on the lattice. A procedure is outlined for 
handling such finite-volume corrections, which must be applied directly to the 
electric form factors themselves rather than to the charge radii. 
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1. Introduction 

Lattice QCD provides important non-perturbative techniques for the analysis 
of many observables. One of the notable features of lattice QCD is that it must 
be performed in a finite volume. The associated effects can be used to access in- 
teresting physics. For example, multi-hadron states are only resolvable at finite 
lattice sizes; the discrete energy eigenvalues becomes increasingly close together 
as the box size becomes large. The finite-volume nature of lattice QCD has in- 
teresting consequences, some of which require careful attention. For example, 
although regularization in both the infrared and ultraviolet regions is an automatic 
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feature of lattice QCD with a finite lattice spacing, finite-sized phenomena, such 
as the virtual pion clouds that surround hadrons, become distorted. This results in 
deviations in the values of lattice observables that can become significant in the 
chiral regime [111, |2|]. Therefore, a method for correcting finite-volume effects by 
estimating their size is sought, using a complementa ry appr oach, such as chiral 
effective field theory (^EFT) [Z,2,A,^^1,3.^M ■ 

Study of the quark mass dependence of lattice QCD simulation results can 
be particularly insightful for examining the chiral properties of hadrons. In order 
to compare lattice simulation results with experiment, an extrapolation to infinite 
volume is additionally required. In order to relate lattice calculations to exper- 
imental results, it is essential to incorporate the low-energy features of QCD in 
order to obtain reliable extrapolations in both quark mass and volume. 

In lattice QCD, form factors are measured at discrete momentum transfer, 
corresponding to the quantization of the momentum modes on the finite spatial 
volume a §,0, 0, Q E [H H 0, S Q . Once form factors have been ex- 
tracted from lattice simulations, they are typically converted directly into charge 
radii. The essential difficulty lies in the definition of the charge 'radius' at finite 
volume (more precisely, the slope of the electric form factor). In order to define 
the radius, a derivative must be applied to the electric form factor, with respect 
to a small momentum transfer, = - ql- This approach breaks down on the 
lattice, where only discrete momentum values are allowed. 

In most cases, calculating the finite-volume corrections poses no essential 
problems [|2. 11. 12, 22, 23, 24, 25]. However, because of the absence of a contin- 
uous derivative on the lattice, the treatment of the electric charge radius is more 
subtle [0,l26l]. Therefore, a method is outlined for handling finite-volume correc- 
tions to a given lattice simulation result. 

Finite- volume charge radii are calculated using the finite-volume electric form 
factors, G|(2^), with an allowed value on the lattice. An extrapolation func- 
tion in will also be introduced to enable a derivative to be taken. It will be 
shown that the finite- volume corrections to the loop integrals must be applied be- 
fore the conversion from the form factor to the charge radius. By applying the 
finite-volume corrections directly to the electric form factor, and ensuring that the 
procedure preserves the electric charge normalization, an extrapolation in is 
then chosen in order to construct an infinite-volume charge radius, defined in the 
usual manner. 
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2. Effective field theory 



In heavy-baryon chiral perturbation theory (^PT), it is usual to define the 
Sachs electromagnetic form factors, Ge,m, which parametrize the matrix element 
for the quark current, 7^, as 

{B(p')\J^\B{p)) = 

u'ip') L GM) + '"'^^^"'^'^'^'^ Gm{Q')\ u\p), (1) 

where v is the velocity and is the positive momentum transfer = -q^ = 
~(p' ~ pf- Lattice QCD results are often constructed from an alternative repre- 
sentation, using the form factors Fi and F2, which are called the Dirac and Pauli 
form factors, respectively. The Sachs form factors are simply linear combinations 
of Fi and F2 

Ge(Q^) = FiiQ^) - ^F2iQ^) , (2) 
Ami 

GMiQ^) = FiiQ^) + F2iQ^). (3) 

In the heavy-baryon formulation of Eq. the spin operator, Sy = -^Jsiy^, y^JVy, 
is introduced HH- 

The external-momentum-dependent electric form factor, Ge{Q^), allows a charge 
radius to be defined in the usual manner 

(r )e = lim -6 ■ (4) 

2.1. Loop integral definitions 

The loop integrals (in the continuum limit) that contribute to the electric form 
factor are invariant under arbitrary translations of the internal momentum k — > 
k + 6k. However, a finite- volume sum over discrete loop momenta is only invari- 
ant if 6k is an allowed value of momentum on the lattice. The loop integrals in 
the heavy-baryon approximation that correspond to the leading-order diagrams in 
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Figs. [T] through [3] are obtained by performing the pole integration for ko 



n J < 



k-4^ 



k k-q 



(5) 

(6) 
(V) 



where co^ = yjk^ - ml, m„ is the pion mass, and A is the Delta-nucleon mass- 
splitting. The chiral coefficients ;\f/v, ;\fA and;^, are derived from couplings arising 
in the Lagrangian of chiral perturbation theory [|29ll 



prot 

Xn ~ 



prot 



xr = 



{D + Ff = -XT\ 



4C^ 



1 



Xa 



= -Xt 



(8) 

(9) 
(10) 



where the value of the pion decay constant is f,^ = 92. 4 MeV. The values for the 
couplings are estimated from the SU(6) flavor-symmetry relations jlSi BOt] and 
from phenomenology: D = 0.76, F = |D and C = -2D. 

2.2. Finite-volume corrections 

The finite- volume correction for a loop integral is defined through the conven- 
tion, following AUton et al. 113 ill 



y - fd^ 



J, 



(11) 



for integrand I. Note that in estimating the size of the finite- volume corrections 
using this method, a lattice spacing parameter is not introduced. 

In a finite-volume calculation of the electric charge radius, the first challenge 
involves the definition of the radius in terms of the derivative in Eq. dH). Since only 
certain, discrete values of momenta are allowed on the lattice, the derivative may 
only be constructed from these allowed momenta when calculating finite-volume 
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Figure 1: The pion loop contributions to the electric charge radius of a nucleon. All charge con- 
serving pion-nucleon transitions are imphcit. 
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Figure 2: The pion loop contribution to the electric charge radius of a nucleon, allowing transitions 
to the nearby and strongly-coupled A baryons. 
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Figure 3: The tadpole contribution at 0{mq) to the electric charge radius of a nucleon. 
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corrections. This crucial observation becomes apparent when a comparison is 
made between the loop integrals evaluated at allowed, and unallowed, values of 
momentum transfer q on the lattice. The comparison is shown in Figs. 4 and 5, 
in which the finite- volume correction to the one-pion loop integral (Eqs. ([5]) & 
©) is plotted as a function of box size L. In Fig. 4, the momentum transfer, 
q, is taken to be ^^in = 2nn/L (where n is a lattice unit vector). A momentum- 
translation of 6k = Kniiii/2 is not an allowed value, and the finite- volume correction 
is inconsistent. In Fig. 5, 6k = k^i^ is an allowed value, and the translated and 
untranslated results for the finite-volume correction are identical. This is simply 
a consequence of spoiling the continuous symmetry by the discretization of the 
momenta. 

To obtain the integrals (T) that contribute directly to the electric charge ra- 
dius, one simply takes the derivative of each infinite-volume integral T{Q^) with 
respect to momentum transfer as — > 0. 

Note that the procedure for calculating the finite-volume corrected electric charge 
radii uses only the infinite- volume integrals (T). Fitting methods are then more 
easily applied at infinite volume. Therefore, the derivative in Eq. (fT2l) need not be 
discretized, but may remain a continuous derivative. 

3. Procedure for obtaining the charge radius at finite volume 

Ideally, one should start with the form factors as extracted from the lattice, and 
only convert them into radii (using a suitable Ansatz for the 2^ -behaviour) after 
correcting for lattice finite- volume effects. At large L, and hence with numerous 
small values, such an Ansatz may be improved upon by a formal expansion 



from;^fEFT, such as that used in Ref. [|32ll . There are many possible choices of 



Ansatz in relating the electric form factor to the charge radius. For example, a 
dipole Ansatz yields the following formula: 

G^(eL) = (i + ^^^^) , (13) 

where = q^ - (E^ - M^)^. Another example is the inverse quadratic with 



two fit parameters, as used by Collins et al. [|21l1 
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Figure 4: (color online). Finite-volume correction for the leading-order loop integrals contributing 
to Ge, with q - ^min. The choice of 8k - k^^ml'^ is not an allowed value on the lattice. The 
momentum translated and untranslated behaviour of the finite-volume correction are inconsistent 
with each other 
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Figure 5: (color online). Finite-volume correction for the leading-order loop integrals contributing 
to Ge, with q - ^min- The choice of 6k - k^^^ is an allowed value on the lattice. Therefore, the 
momentum translated and untranslated behaviour of the finite-volume correction is identical. 
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with the charge radius obtained through {r )e = 6 a. 

One needs to transform the finite- volume form factors, Gp(Q^ ), to infinite- 
volume before the charge radii are calculated. The finite-volume correction to 

is achieved by subtracting the difference between the sum and integral loop 
contributions, from Eq. (fTTI) 

G~(e') = G^Ce^) - 6dT(Q^)], (15) 

where T = Tn + 7~a + Ttad- In order to maintain conservation of charge, a subtrac- 
tion of the finite- volume corrections at zero momentum transfer, 6[T{Q^ = 0)], is 
introduced 

SdTiQ^)] ^ SdTiQ^)] - SdnO)]. (16) 

The second term of Eq. (fT6l) ensures that both infinite- and finite- volume form fac- 
tors are correctly normalized, i.e. G^(0) = 1 = G^(0), independent of the nucleon 
momentum. While this normalization procedure assumes that charge conserva- 
tion is satisfied in a finite volume, the following numerical and ;tEFT analyses 
confirm that this is indeed so. 

A lattice QCD calculation was undertaken to test charge conservation in the 
rest frame, and in several boosted frames of increasing momenta. Using a con- 
served vector current and the prescription outlined in Ref. [il6i1 . the Sachs electric 
form factor G^ at zero momentum transfer was calculated for p = Q, k^m, 2 ^min 
and 3 fc^in, where k^^^ is the minimum available three-momentum on the lattice. In 
all cases considered, the extracted value of G^{Q^ = 0) was consistent with unity 
up to the level of precision present in the propagators used in the calculation, 
indicating that charge conservation is satisfied. 

To provide support for this calculation from an effective field theory stand- 
point, the electromagnetic Ward Identity was verified in;^fEFT with a renormalized 



current, using the same technique as Djukanovic et al. [|33|l . Since the tree-level 
nucleon-photon vertex r^^tree ~ ^T^T'^ automatically satisfies the Ward Identity, 
it is sufficient to show that 

ip' - P)Xn = [^Nif) - ^Nif)] , (17) 

where p' and p are the outgoing and incoming momenta, respectively. I^n de- 
notes the standard one-pion loop contributing to the mass of the nucleon. 
denotes the one-pion loop where all possibilities of a one-photon attachment are 
summed. For example. Fig. [U shows the contribution to the electromagnetic Ward 
Identity for a photon-pion vertex. The diagrams with photon attachments to the 
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internal/external baryon, and at the baryon-pion vertex, must also be included. In 
proving the Ward Identity, one does not make use of any infinite-volume proper- 
ties of the loop integrals. One may freely replace the three-dimensional part of 
each integral with a sum, using the prescription in Eq. (fTTl) . without disrupting the 
conservation of charge. Thus, the electric charge is conserved at finite volume in 
lattice QCD, and to one-loop order in;tfEFT. Of course, the effective field theory 
will continue to respect this symmetry to all orders. 

The infinite-volume charge radii (r^)^ can now be recovered from the form 
factors by using an extrapolation in Q^. For example, in the case of the dipole 



Ansatz 



. 12 



1 



(18) 



Once the infinite-volume charge radii have been obtained, one then has best es- 
timates for charge radii at any given set of lattice parameters. Importantly, this 
provides a value of the charge radius at different quark (pion) masses, and hence 
enables the use of continuum ^EFT to fit the quark mass dependence. Details of 
such an extrapolation will appear in a forthcoming paper. 

As a demonstration of the method, lattice QCD results from QCDSF Collabo- 
ration [21] for the isovector nucleon electric charge radius are corrected to infinite 
volume, and shown in Fig. [3l The lattice calculation uses the two-flavor 0(a)- 
improved Wilson quark action. The points displayed satisfy m„L > 3, and the box 
sizes for each point range from 1.7 to 2.9 fm. The infinite volume form factors are 
calculated via 

G^(e') = G^(e') - (sdTiQ^)] - sdrm), m 

with TiQ^) = + TA(e') + TUQ^) (20) 



In the construction of the charge radii, the ^^-extrapolation Ansatz of Eq. (fT4)) is 
used, for example. Evidently, the infinite-volume points are larger in radius than 
the finite-volume points, and closer to the experimental value of (r^Y^"" = 0.88 



fm^ [34L135] 



4. Conclusion 

Direct finite-volume corrections to charge radii are not well defined on the 
lattice. The use of continuous derivatives in constructing the electric charge ra- 
dius leads to inconsistent results in the finite-volume corrections. Alternatively, 
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Figure 6: (color online). Lattice QCD results for the isovector nucleon electric charge radius from 
QCDSF Collaboration at their original box sizes, and corrected to infinite volume. 

a satisfactory definition of radii can be achieved by implementing finite-volume 
corrections directly to the electric form factors, as evaluated at discrete values of 
q^. Subsequently, the resultant finite- volume-corrected form factor may then be 
converted into a charge radius using an appropriate extrapolation in Q^. This then 
defines a charge radius suitable for comparison with experiment and continuum 
^EFT analysis. 

This research is supported by the Australian Research Council through grant 
DP110101265. 
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